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I HTRODUCT I  OH 


This  Final  Report  suaaarlzes  the  principal  results  of  research  activities 
undertaken  under  Rray  Research  Office  Grant/Contract  Ho.  DRRLQ3-89-K-0134,  1  Jut 
1389  -  30  Septeaber  1990,  under  the  auspices  of  the  Research  Foundation  of  the 


State  University  of  Hee  York.  Additional,  tore  detailed  eater ial  is  included  in 
the  aanuscripts  cited  in  the  Bibliography  as  entries  [1]-[8],  each  derived  in 
part  froa  the  support  and  activities  of  this  contract. 

The  Principal  Investigator  Peter  0.  Scott,  Departaent  of  Electrical  and 
Coaputer  Engineering  (ECE),  State  University  of  Nei  York  at  Buffalo  (SUHVB) ,  «as 
aided  in  this  project  by  several  able  individuals.  Mr .  Li  fan  Hua,  candidate  for 
the  Ph.D.  degree  in  ECE  SUHYB,  assisted  in  ail  aspects  of  the  aork  and  aas 
supported  under  the  contract  for  a  brief  period  of  tlae.  It  is  anticipated  that 
he  sill  coaplete  requireaents  for  the  Ph.D.  degree  within  the  next  six  aonths, 
partly  as  a  result  of  his  activities  and  support  under  this  contract.  Hr.  Gongeei 
Xie,  ECE  SUHYB  Research  Rs*nciate,  contributed  extensive  efforts  at  the  optical 
bench  aithout  any  charge  to  the  contract.  Dr.  David  T.  Shae,  Professor  of  ECE 
SUHYB,  offered  his  aerosol  laboratory  facilities  and  supplies,  along  eith 
frequent  technical  consultations  and  suggestions,  also  aithout  any  charge  to  the 
contract.  Uithout  the  aagnaniaous  support  of  the  these  tao  individuals,  scope  and 
activities  sou  Id  have  been  auch  attenuated. 

As  aore  coaplete I y  described  in  RRO  Proposal  Huaber  27078-GS,  the  central 
focus  of  the  research  activities  supported  under  this  contract  aas  investigation 
of  the  feasibility  of  utilizing  Fresnel  holography,  together  eith  digital 
reconstruction  techniques,  to  assess  the  pair  correlation  of  aodel  particle 
systeas.  Experiaental  activity  consisted  in  the  fabrication  of  an  optical  cell  in 
ehich  the  particles  could  be  suspended  and  of  the  particle  systeas  theaselves, 
alignaent  of  the  YAG  laser  and  optical  bench  for  inline  Fresnel  holography,  and 
acquisition  and  processing  of  the  holograas  of  the  optical  cell  contents. 
Representative  sections  of  the  processed  holograas  aere  subsequently  digitized 
under  high  aagni f icat ion,  input  to  a  general  purpose  coaputer,  and  analyzed  using 
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algorithms  me  developed  to  estimate  each  visible  particle's  centroid  location  at 
the  highest  attainable  resolution.  The  resulting  three  dimensional  data  base  of 
particle  location  parameters  mas  mapped  into  an  pair  correlation  function 
estimate  using  standard  non-parametric  distribution  function  estimation 
technique. 


BACKGROUND 

By  identifying  interparticle  separations  of  higher  and  lomer  probability  than  the 
completely  random  case,  the  pair  correlation  function  (or,  In  the  isotropic  case 
the  radial  distribution  function)  reveals  preferred  spacing  intervals  even  in 
materials  such  as  gases  and  fluids  that  have  no  obvious  long  range  order.  These 
preferred  spatial  scales  predict  mavelength  dependent  scattering  and  attenuation 
effects  for  move  structured  energy  propagating  through  the  particle  field. 
Anisotropic  pair  correlation  suggests  dependence  on  the  direction  of  move 
propagation  relative  to  the  orientation  distribution  of  the  particle  ensemble  as 
•ell  as  mavelength  dependence. 

The  pair  correlation  is  a  primary  input  parameter  in  recent  predictive 
multiple  scattering  theories  of  move  propogation  through  aerosols  and  hydrosols. 
In  the  Uaradan  approach  [9]-[11],  the  quasicrysal I ine  approximation  is  used  to 
characterize  ensemble  field  coefficients  in  terms  of  the  pair  correlation,  from 
•hich  a  dispersion  equation  is  derived  using  Tmersky's  extinction  theorem.  This 
method  for  prediction  cf  scattering  and  attenuation  of  »ave  energy  is  of  broad 
general ity,  appluing  to  Lulh  electromagnet i c  and  compression  eaves  in  air,  mater, 
vacuum  and  other  media. 

Techniques  employed  to  estimate  the  pair  correlation  of  "interesting" 
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particle  fields  include  truncations  based  on  theoretical  arguments  (such  as 
Perkus-Veoich  and  hypernetted  chains),  digital  computer  s leu  I  at  ion  employing 
Monte  Carlo  or  aolecular  dynamics,  and  measurement  of  x-ray  structure  factor. 
Truncations  of  the  equations  of  state  yield  systems  mhich  can  only  be  soloed,  or 
eoen  reliably  esti sated,  in  the  simplest  cases  such  as  monodisperse  dilute  hard 
spheres.  Lagos  and  Seville  [12]  have  recently  soloed  the  Perkus  -  Veolch  equation 
numerically  for  hard  spherocyl inders  (particles  formed  by  capping  right  circular 
cylinders  mith  half-spheres).  They  failed  to  find  eoen  rough  corroboration  for 
current  pair  correlation  mixture  estimates  based  on  perturbation  theories  or 
simplified  integrals. 

Digital  simulations  are  much  more  general  in  scope  and  model  parameters, 
but  for  realistic  polydisperse  dense  particle  fields,  simplifying  assumptions 
must  be  applied  mhich  bring  the  results  into  question,  and  eoen  then  computing 
time  becomes  a  cost  problem  eoen  on  modern  high  speed  digital  computers.  This  is 
particularly  evident  in  the  case  of  relatively  high  density  non-spherical 
particles  [13].  X-ray  structure  factor  measurement  Is  an  attractive,  direct 
method  for  pair  correlation  estimation,  but  is  limited  to  cases  mhere  particles 
are  of  the  same  granularity  as  the  x-ray  save  length.  Thus  the  need  exits  for  ' 
method  of  assessment  of  pair  correlation  mhich  mould  complement  the  capabilities 
of  the  existing  techniques.  This  is  the  opportunity  addressed  by  the  holographic 
method  for  pair  correlation  measurement  in  model  systems  investigated  in  the 
present  research  project. 


PROBLEM  STATEMENT 

The  central  problem  to  be  addressed  is  the  accurate  assessment  of  the 
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distribution  of  centroids  of  every  particle  • i t h i n  a  characteristic  test  volute 
of  the  optical  cell.  Underlying  this  tensuration  goal  is  the  hypothesis  that  the 
ensetble  of  todel  particles,  scaled  in  size  relative  to  the  target  particle  field 
•hose  pair  correlation  is  to  be  esti sated,  till  accurately  reproduce  the  (scaled) 
pair  distribution  function  of  the  target  paticie  field  those  pair  correlation  is 
being  estimated.  Note  that  this  assumption  is  plausible  only  in  the 
non-interacting  hard  particles  case  (in  which  the  only  interparticle  interaction 
consists  of  a  non-interpenetrability  constraint).  Other  particle-particle 
interactions,  such  as  electrostatic  or  magnetic  coupling,  would  be  tore  difficult 
to  situ  I  ate  in  the  todel  particle  field  with  accuracy,  since  in  general  such 
forces  do  not  scale  coherently  with  particle  dimension. 

The  three  dimensional  particle  centroid  location  database  must  be  acquired 
•ithin  a  fet  seconds  after  the  todel  particle  field  is  last  stirred,  in  order  to 
prevent  settling  and  distortion  of  the  particle  separation  and  orientation 
geoietry.  Each  particle  in  a  sample  portion  of  the  optical  cell  tust  be  observed 
and  its  centroid  computed  accurately.  The  accuracy  of  centroid  location,  i.e.  the 
resolution  of  the  measurement  system,  must  be  significantly  greater  than  a 
standard  particle  length  in  order  to  detect  and  locate  narrow  features  (peaks  and 
nulls)  in  the  pair  correlation  function  with  useful  accuracy.  Finally,  the 
resulting  list  of  particle  locations  in  three  space  must  be  mapped  into  a  pair 
correlation  estimate,  taking  into  account  the  resolution  -  smoothness  tradeoff 
inherent  in  all  probability  density  function  estimation  methods.  The  estimation 
method  must  have  sufficient  resolution  to  detect  significant  fluctuations  from 
the  purely  random  distribution  case  even  when  they  are  relatively  narrow,  while 
preserving  as  low  an  estimate  variance  as  is  attainable  consistent  with  the 
required  resolution. 
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METHODS 


The  optical  cell  used  mas  approx  I lately  2  ci  on  a  side,  lade  of  high 
quality  optical  glass,  iith  opposing  surfaces  carefully  aligned  and  glued  for 
para  I  lei inity.  The  particles  vere  standard  latex  spheres  of  various  sizes  ranging 
froi  5  eicroieters  to  100  micrometers,  as  sell  as  cylindrical  iron  fibers 
produced  in  our  laboratory,  licrotoied  and  sieved  to  size.  The  eedia  in  ihich  the 
particles  eere  suspended  eas  high  purity  later,  fill  experiments  mere  done  at  rooi 
temperature. 

The  optical  setups  are  shorn  in  Figures  1  and  2.  Figure  1  shoes  the  in-line 
hologram  acquisition  configuration,  in  ehich  the  optical  elements  are  i 1 1 um i anted 
by  a  VflG  laser,  fl  high  poeer  pulsed  laser  eas  utilized  to  permit  shortening 
exposure  time  sufficiently  so  that  high  resolution,  relatively  Insensitive 
hologram  emulsions  could  be  used  and  yet  avoid  motion  blur.  The  optical  cell  mas 
located  at  position  D.  Figure  2  shoes  the  monitoring  configuration  used  to  check 
digital  reconstructions  by  direct  optical  measurement.  The  conventionally 
reconstructed  image  is  routed  to  to  TU  monitor  and  via  analog  to  digital 
conversion  (ROC)  to  a  general  purpose  computer  for  comparison  with  digital 
reconstruction  (DEC  URX  780  and  SUN  3/260  computers  mere  used).  The  monitor  mas 
useful  in  surveying  the  holograms,  tuning  the  digital  algorithms  and  verifying 
their  basic  particle  counts  along  mith  coarse  verification  of  the  centroid 
computations. 

The  digital  algorithms  used  to  analyze  the  contents  of  the  R0  converted 
holograms  are  based  on  the  tmin-image  elimination  methods  previously  developed  in 
this  laboratory  [ 14]— [ 15] .  In  order  to  maximize  the  resolution  (thus  reducing  the 
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possibility  of  aissing  closely  opposed  particles  and  reducing  the  centroid 
location  error  variance  of  the  enseable)  a  "superresolution"  signal  processing 
layer  was  built  upon  this  base  [1],  [2].  In  this  technique  the  bandeidth 
limitation  of  the  reconstruct  ion  effectively  iaposed  by  the  physical  size 
(holograa  plate  foraat)  limitation  is  stretched  by  an  iterative  spectral 
extrapolation  procedure.  The  effect  is  to  computationally  recreate  a  portion  of 
the  hologram  that  mas  lost  at  the  borders  of  the  holographic  plate  thereby 
permitting  higher  frequency  components  (i.e.  higher  resolution)  in  the  object 
field  reconstruction.  This  algorithm  is  designated  as  the  Phase  Retrieval 
Spectral  Continuation  (PRSC)  algorithm.  Final ly,  the  mapping  from  detected 
centroid  locations  in  three  space  to  a  final  pair  correlation  function  estimate 
mas  made  using  the  standard  Parzen  function  method. 

RESULTS 

R  central  result  of  this  study  mas  the  successful  upgrading  of  algorithms 
originally  designed  for  general  particle  field  inline  hologram  decoding  to  the 
present  purposes.  The  resolution  of  the  algorithms  presented  in  [M]  and  [15] 
needed  to  be  improved  to  permit  sufficiently  accurate  and  consistent 
determination  of  the  centroid  locations  of  each  of  the  particles  in  the 
holographic  test  region. 

The  referenced  algorithms  effectively  reproduced  by  digital  means  the 
resolution  available  mith  high  quality  conventional  (optical)  resolution  once  the 
tmin  image  artifact  had  been  suppressed.  This  required  enhancement  for  effective 
pair  correlation  estimation  for  tmo  reasons.  First,  the  reduced  axial  resolution 
inherent  in  in- line  holography  [16]  together  mith  the  shadoming  of  particles 
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located  in  a  narrow  cone  behind  tore  proxiaal  particles  (relative  to  boresight 
distance  froa  the  laser  source)  aade  it  possible  that  deeply  shadowed  particles 
aight  not  be  detected.  This  effect  aust  be  ainlaized,  since  it  leads  to 
undercounting  of  closely  opposed  particles  which  happen  to  be  aligned  along  the 
optical  axis,  and  thus  to  systeaatic  bias  reducing  the  pair  correlation  at  low 
separations.  Second,  the  error  uariance  of  the  pair  correlation  estiaate  at  given 
separation  is  directly  proport ionai ,  in  a  given  run,  to  the  individual  variance 
in  the  estiaate  of  centroid-centroid  separation  for  a  given  particle  pair  falling 
within  that  aeasuraent  bin  reduced  by  the  square  root  of  the  nuaber  of  particle 
pair  separations  within  that  bln.  Since  there  are  relatively  few  particle  pairs 
at  the  saaller  separations,  the  variance  is  highest  exactly  where  the  pair 
correlation  is  of  aost  interest  and  aust  be  reduced  as  auch  as  possible  in  order 
that  interesting  features  not  be  lost  under  the  aeasuraent  noise. 

This  enhanceaent  was  achieved  by  using  superresolution  signal  processing 
aethods  described  in  [1]  and  [2].  fln  increase  of  factors  of  two  to  eight  was 
aeasured  in  the  available  resolution,  dependent  on  experiaental  signal  to  noise 
ratio.  To  deteraine  the  likelihood  of  aissed  particles,  sieulations  were 
eaployed.  It  was  found  that  only  under  very  low  probability  circuastances,  such 
as  the  coaplete  shadowing  of  a  saaller  particle  within  one  to  two  interparticle 
separations  directly  behind  the  larger  one  was  a  particle  likely  to  be  aissed. 
Thus  at  all  particle  densities  consistent  with  accurate  in-line  aeasureaent,  i.e. 
at  least  75X  of  the  illuaination  reaching  the  holograw  along  the  undiffracted 
reference  pathway  [16],  the  resolution  of  the  resulting  technique  was 
satisfactory  to  the  present  purpose.  Studies  of  both  latex  spheres  and 
cylindrical  particle  populations  at  low  densities  (up  to  0.05X  voluae  fraction 
density)  indeed  verified  this  assertion.  Uithin  experiaental  error  the  resulting 
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pair  correlation  estimates  had  the  familiar  "veil  stirred"  configurations. 

Unfortunately,  the  optical  cells  aual table  for  this  project  proved 
unsuitable  for  higher  density  trials.  In  order  to  preserve  the  751  criterion,  a 
straightfoaard  calculation  shoes  that  at  the  volume  fraction  of  IX,  the  border  of 
interestingly  "turgid"  fields,  it  is  necessary  that  the  cell  be  only  on  the  order 
of  ten  particle  lengths  thick  (measured  along  the  optical  axis;  the  other 
dimensions  may  be  arbitrarily  large).  The  narrowest  cell  available  was  2  ca 
thick,  requiring  particles  too  large  to  be  effectively  holograaaed.  Our  efforts 
to  fabricate  cells  aith  thickness  1-2  aa  aere  unsuccessful  aithln  the  tiae  and 
cost  constraints  (the  fabrication  effort  eas  unsupported  in  this  contract). 
Difficulties  included  the  very  narroa  separation  distances  and  the  high  degree  of 
para  I  lei inity  required  to  prevent  coherent  reflection  artifact. 

CPflCUJSIPHS 

The  principal  conclusion  of  this  feasibility  study  is  that  in-line 
holography,  supported  by  digital  reconstruction  and  resolution  enhancement,  Is  a 
viable  complement  to  analytical  approximations,  digital  simulations  and  x-ray 
structure  factor  measurement  for  the  determination  of  pair  correlation  aithin 
some  model  aerosol  systems,  fit  lorn  densities  the  me  1 1  stirred  solution  is 
recovered  aithin  experimental  error,  fit  higher  densities  very  narroa  optical 
cells  are  required  to  maintain  the  quality  of  the  in-line  holograms.  While  no 
experimental  data  mas  available,  all-digital  simulations  suggest  the 
superresolution  algorithms  should  be  equally  appicable  in  the  high  density  as  in 
the  lorn  density  regime. 


9 


B I BL I OGHRPHV 

[1]  L.  Hua,  0.  Sham,  P.D.  Scott,  "fin  adaptive  constraint  algorithm  for  image 
understanding  of  digital  reconstructed  in-line  holograms,"  Proceedings  of  the 
1990  IEEE  International  Conference  on  Systems,  Han  and  Cybernetics,  pp. 

150-6. 

[2]  L.  Hua,  G.  Xie,  0.  Shai,  P.0.  Scott,  “Resolution  enhancement  in  digital 
in-line  holography,'  Proceedings  of  the  SPIE,  accepted  for  publication  to 
appear  1991. 

[3]  G.U.  Xie,  L.F.  Hua,  S.  Patel,  P.  Scott,  D.T.  Sham,  "Double  exposure 
holography  for  fast  moving  fibers  mith  submicron  diameter,*  Rpplied  Optics, 
accepted  for  publication  to  appear  1991. 

[4]  Gongmei  Xie,  Peter  Scott,  David  T.  Sham,  Vimo  Zhang,  “Particle  measurement 
using  overlapping  optical  field  technique,"  under  reviem  by  Rpplied  Optics. 

[5]  Gong-mei  Xie,  Peter  Scott,  David  T.  Sham,  Vimo  Zhang,  "R  nem  optical 
technique  for  particle  sizing  and  velocimetry, “  Optics  Letters,  accpeted  for 
publication  to  appear  1991. 

[6]  G.U.  Xie,  L.F.  Hua,  S.  Patel,  P.  Scott,  D.T.  Sham,  "Using  dynamic  holography 
for  iron  fibers  mith  submicron  diameter  and  high  velocity,"  Proceedings  of 
the  SPIE,  accepted  for  publication  to  appear  1991. 

[7]  G.U.  Xie,  L.F.  Hua,  S.  Patel,  P.  Scott,  D.T.  Sham,  "Single  shot  holography  of 
fast  moving  fibers  mith  subm.cron  diameter,"  under  reviem  by  Optical 
Engineering. 

[8]  G.  Xie,  P.  Scott,  "Interpreting  digital  nolograms,"  Proceedingo  of  thi  1991 
International  Congress  on  Optical  Scier-ce  and  Engineering,  The  Hague, 
Hetherlands.  Accepted  for  publication,  to  appear  1991, 

[9]  U.U.  Uaradan,  U.H.  Bringi,  U.K.  Uaradan,  "Coherent  electromagnetic  move 


10 


propogation  through  randomly  distributed  dielectric  scatterers,"  Physics 
Reviem  D  vol.  19  (1979)  p.  2480. 

[10] U.U.  Uaradan,  V.  11a,  U.K.  Uaradan,  "Propogator  todel  including  multipole 
fields  for  discrete  random  media,"  Journal  of  the  Optical  Society  of  fleer ica 
A  vol.  2  (1985)  p.  2195. 

[11] U.U.  Uaradan,  U.K.  Uaradan,  V.  da,  "Effects  of  nonspherical  statistics  on 
electromagnet ic  move  propogation  in  discrete  random  media,"  Journal  of  the 
Optical  Society  of  fleer ica  fl,  uo I .  2  (1985)  p.  491. 

[123S.  Lagos,  P.  Sevilla,  "Solution  of  the  Perkus  Vevich  equation  for  hard 

spherocy I inders  I.  The  entire  pair  correlation  function,"  Journal  of  Chemical 
Physics,  vol.  89  (1988)  p.  4349. 

[13] 0.  Rigby,  fl.  floe,  "Molecular  dynamics  simulation  of  polymer  liquid  and  glass 
II.  Short  range  order  and  orientation  correlation,"  Journal  of  Chemical 
Physics,  vol.  89  (1988)  p.  5280. 

[14] L.  Onural,  P.  Scott,  "Digital  decoding  of  in-line  holograms,"  Optical 
Engineering,  vol.  34  (1987)  p.  2287. 

[15] li.  Liu,  P.  Scott,  "Phase  retrieval  and  to  in- image  elimination  for  in-line 
Fresnel  holograms,"  Journal  of  the  Optical  Society  of  America  fl,  vol.  4 
(1987)  p.  1118. 

[16] P.  Hariharan,  Optical  holography.  Cambridge  University  Press,  London,  UK, 
1984. 


YAG  Laser 


Ml.  M2.  M3.  M4.  M5.  M6: 
Mirrors  far  high  power  laser 

SI.  S2:  Splitter 

1:  bis  _ 


LI:  Negative  lens 
L2.  L3;  line  • 

D:  Disperser 


L4:  Microscope  lens.  25x.  0.31  NA 
( long  working  distance) 

HP:  Holographic  plate 


_ _  AS:  Allgnnlng  System  Ma: 

Figure  1.  Holograa  acquisition  optical  bench  setup. 


He-Ne  Laser 


Camera 

13  *  12 

U 

z  *< 

TV  Modter 

_ f=a 

LI  :  Microscope  Lens 

PH:  Pin  Hole 

L2  :  Collimator  Lens 

A  Attenuator 

1111 

L3  :  Mag  Hying  Lens 

H  :  Hologram 

llll  — 

Digital  Image 
Processing  System 

Ml.  M2  :  Mirror 

Figure  2.  Optical 

reconstruction  and  aonitorlng. 

RPPEMD1X 


The  first  pages  of  aanuscripts  [1]-[8]  listed  in  the  Bibliography  follow  in 
order.  Each  of  these  aanuscripts  «as  based  in  part  on  the  activities  and 
support  of  flRO  Grant/Contract  DRRLQ3-89-K-0134.  The  brief  duration  of  the 
contract  (12  aonths  with  a  three  aonth  extension)  explains  why  aost  of  these 
aanuscripts  are  under  review,  or  accepted  but  not  yet  appeared,  at  the  tiae 
of  closure  of  this  Final  Report. 


13 


An  adaptive  constraint  algorithm  for  image  understanding 
of  digitally  reconstructed  in-line  holograms 
by 

Lifan  Hua,  David  T.  Shaw  and  Peter  D.  Scott 
Dept,  of  Electrical  and  Computer  Engineering 
State  University  of  New  York  at  Buffalo 
Bell  Hall,  Amherst  NY  14260 

Abstract 

Accurate  quantitative  decoding  of  the  3D  information  encoded  in  in-line  holograms,  at  best  difficult 
optically,  may  be  achieved  by  A/D  sampling  of  the  hologram  followed  by  digital  reconstruction.  Using  standard 
reconstruction  technique,  the  resulting  image  bandwidth  is  limited  by  the  size  of  the  hologram.  In  addition,  the  phase 
ambiguity  inherent  in  magnitude-only  hologram  recording  yields,  along  with  the  desired  object  reconstruction,  an 
out-of-focus  conjugate  artifact  called  the  twin  image.  Both  limited  bandwidth  and  twin  image  restrict  the  available 
resolution  well  below  the  theoretical  diffraction  limit  In  this  paper  an  algorithm  is  presented  which  addresses  both 
problems  by  iteratively  combining  phase  retrieval  and  spectrum  continuation  to  produce  estimates  of  (he  phase  of  the 
recorded  hologram,  and  both  magnitude  and  phase  of  the  hologram  beyond  its  physically  recorded  boundaries. 

Since  algorithms  based  on  spectral  continuation  are  sensitive  to  the  constraints  imposed  on  the  extent  of 
the  objects  being  imaged  (which  must  be  space  limited),  a  method  for  selecting  these  constraints  adaptively  was 
developed.  Effective  use  of  adaptive  constraints  greatly  accelerates  convergence  of  the  algorithm  compared  to  fixed 
constraints.  The  chosen  adaptive  constraint  selection  rule  is  shown  to  have  an  error-correcting  property  in  which 
over-constraint  is  naturally  corrected  and  the  constraint  set  boundaries  stably  relax  toward  the  actual  object 
boundaries.  Examples  demonstrate  the  degree  of  improvement  in  resolution  over  conventional  reconstruction,  and 
over  reconstruction  with  phase  retrieval  only.  This  algorithm  suggests  the  possibility  of  super-resolution 
holography,  in  which  the  diffraction  limit  is  exceeded.  While  difficult  to  achieve  in  practice,  super-resolution  has 
been  demonstrated  for  other  kinds  of  imaging  but  not  heretofore  considered  for  holography. 

1.  Introduction 

Under  a  variety  of  conditions  a  multidimensional  signal  may  be  reconstructed  from  partial  information 
about  its  transform,  for  instance  its  Fourier  or  Fresnel  transform!  1-4].  While  the  ability  to  reconstruct  a  signal  from 
magnitude  information  only  would  be  useful  in  several  applications  settings,[5-9]  magnitude  alone  is  not  in  general 
sufficient  to  reconstruct  a  signal.  For  example,  holography  is  a  useful  modality  for  capturing  an  object  set's  three 
dimensional  distribution.  But  only  the  intensity  of  the  optical  field  is  recorded  on  the  film  emulsion  or  other 
photosensitive  medium. 

The  clarity  and  resolution  of  conventionally  reconstructed  holographic  images  is  limited  by  factors 
including  size  of  the  hologram  plate,  grain  and  the  nonlinearity  of  the  holographic  film  and  the  twin  image  artifacts. 
Techniques  for  suppressing  twin  image  in  the  digitally  reconstructed  in-line  holography  by  phase  retrieval  and  by 
modified  inverse  filtering  have  been  reported  [5,6].  These  techniques  do  not  seek  to  overcome  the  reconstruction 
bandlimitation  inherent  in  the  finite  hologram  size.  In  the  case  where  missing  high  frequency  components  is  the 
dominant  resolution  restriction,  the  resolution  of  the  reconstructed  images  will  not  be  improved  much  no  matter  how 
precise  the  estimated  phase  information  is.  In  this  paper,  a  new  method  for  "continuing"  or  extrapolating  a  given 
hologram  is  described.  This  method  employs  a  convolution/deconvolution  iterative  continuation  cycle,  combined 
with  interior  phase  retrieval  since  the  mssing  phase  information  is  required  for  successful  continuation.  This 
algorithm  uses  a  Gerchberg-Saxton  like  "error  energy"  reduction  principle. 

With  standard  constraint  specifications,  considerable  computational  time  is  needed  due  to  the  slow 
convergence  of  the  combined  phase  retrieval  and  spectrum  extrapolation  (phase  retrieval  and  spectrum  continuation 
are  each  known  to  converge  slowly;  here  both  are  being  solved  simultaneously).  By  adaptive  selection  of  the 
constraints  in  the  object  lomain,  the  convergence  may  be  greatly  accelerated  without  inducing  instabilities  or  errors. 

2.  Background  and  problem  formulation 

It  is  assumed  in  this  paper  that  the  objects  are  opaque  and  the  background  is  transparent,  and  the  thickness 
of  each  individual  object  can  be  ignored.  Thus  the  objects  being  imaged  can  be  considered  as  a  set  of  (cross-sectional) 
planar  distribution  function  in  the  3D  object  space.  These  assumptions  are  motivated  by  the  study  of  particle  fields 


14 


Resolution  Enhancement  in  Digital  In-line  Holography 


Lifan  Hua,  G.  Xie,  David  Shaw  and  Peter  Scott 
Department  of  Electrical  and  Computer  Engineering 
State  University  of  New  York  at  Buffalo 
Buffalo,  NY  14260 


Abstract: 

The  image  volume  captured  in  a  in-line  hologram  may  be  reconstructed  optically  or,  following 
analog-to  digital  conversion,  algorithmically  with  the  aid  of  a  digital  computer.  Factors  limiting  the 
clarity  and  resolution  of  the  resulting  reconstructions  (both  optical  and  digital)  include  size  of  the 
hologram,  grain  and  nonlinearity  of  the  film,  and  twin  image  artifacts.  Applications  emphasizing 
display  and  visualization  require  optical  reconstruction,  while  technical  applications  demanding 
detailed  quantitative  image  volume  descriptions  favor  digital  technique. 

An  added  advantage  of  digital  reconstruction  is  the  ease  with  which  postprocessing  may  be 
applied  to  suppress  noise  and  artifact,  or  enhance  resolution.  In  this  communication  a  constrained 
iterative  “super-resolution”  type  algorithm  is  described  and  applied  to  enhance  resolution  of  in-line 
holograms  of  the  particle  field  type.  This  algorithm  is  based  on  the  unique  relationship  between  a 
space-limited  object  and  any  finite  segment  of  its  spectrum.  If  the  truncated  complex  spectrum  of  a 
space-limited  object  is  measured  or  derived,  the  full  spectrum  may  be  continued  in  only  one  way. 

For  in-line  holograms,  only  the  magnitude  of  the  holographic  field  is  recorded.  To  achieve 
super-resolution,  both  magnitude  and  phase  information  are  needed  for  spectrum  extrapolation. 
Therefore,  a  phase  retrieval  scheme  is  joined  with  spectrum  extrapolation  in  the  present  algorithm. 
This  method  uses  a  Gerchberg-Saxton  like  convolution/deconvolution  iterative  cycle  to  trim  the 
iteratively  estimated  data  subject  to  the  constraints  of  the  object  domain  and  hologram  domain.  The 
constraint  in  the  object  domain  is  that  the  object  has  finite  size  in  space,  and  in  the  hologram  domain 
the  constraint  is  the  known  measured  magnitude. 

The  super-resolution  algorithm  has  been  tested  on  both  computer  generated  holograms  and 
real  optical  holograms.  Results  demonstrate  the  ability  to  significantly  improve  resolution  while 
suppressing  the  twin  image,  the  algorithm  is  shown  to  be  robust  against  noisy  data.  The  algorithm 
permits  accurate  and  high  resolved  quantification  of  dynamic  3D  particle  fields. 
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DOUBLE  EXPOSURE  HOLOGRAPHY  FOR  FAST  MOVING 
FIBERS  WITH  SUBMICRON  DIAMETER 

G.W.  Xie,  L.F.  Hua,  5.  Patel,  P.  Scott,  D.T.  Shaw 

Department  of  Electrical  and  Computer  Engineering 
State  University  of  New  York  at  Buffalo,  Buffalo,  New  York  14260 


ABSTRACT 

In  this  paper,  an  optical  technique  which  use  divergent  beam,  high 
magnification  microscope  objective  lens  and  double  exposure  method  to 
record  holograms  of  fast  moving  fibers  (125meters/second)  with 
submicron  diameters  (0.5pm)  is  reported  for  the  first  time.  The  principle 
of  using  divergent  laser  beam  to  protect  high  magnification  microscope 
objective  lens  based  on  ray-tracing  equations  and  the  calculation  of 
magnification  of  whole  process  for  recording  and  reconstructing 
holograms  based  on  optical  theory  are  presented.  An  in-line  double 
exposure  holographic  experimental  system,  utilizing  the  principle  of 
dividing  one  beam  to  two  parts  with  different  optical  distance,  was  set  up 
to  record  the  holograms  of  fast  moving  fibers.  A  set  of  optical 
reconstructions  of  the  holograms  are  given  and  the  measuring  results  are 
presented.  The  results  show  that  by  using  the  method  above  mentioned  , 
the  accurately  measuring  results  of  both  fiber's  submicron  diameter  and 
high  velocity  can  be  determined  at  the  same  time. 
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Particle  Measurement  Using  Overlapping  Optical  Meld  Technique 


Gongwei  Xie,  Peter  Scott,  David  T.  Shaw 

Department  of  Electrical  and  Computer  Engineering 
State  University  of  New  York  at  Buffalo 
Buffalo,  NY  14260 


Yimo  Zhang 


Department  of  Precision  Instrument  Engineering 
Tianjin  University 
Tianjin,  P.  R.  China 


Abstract 

An  optical  scattering  system  is  described  which  uses  the  Overlapping 
Optical  Field  technique  (OOF)  for  simultaneous  measurement  of  particle  size 
and  velocity.  Experimental  results  using  this  system  to  study  rapidly  moving 
small  particles  of  varying  size,  shape  and  density  are  presented. 

Key  Words:  overlapped  optical  field,  light  scattering,  particle 
measurement,  polarized  direction,  measurement  volume, 
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A  New  Optical  Technique  For  Particle  Sizing  And  Velocimetry 

by 

Gong-wei  Xie,  Peter  Scott,  David  T.  Shaw 

Department  of  Electrical  and  Computer  Engineering,  State  University  of 
New  York  at  Buffalo,  Buffalo,  NY  14260 

Yimo  Zhang 

Department  of  Precision  Instrument  Engineering,  Tianjin  University 

Tianjin,  P.  R.  China 

Abstract 


A  novel  Overlapping  Optical  Field  (  OOF)  technique  that  transforms  an 
optical  fieid  of  Gaussian  intensity  into  one  with  uniform  intensity  is 
described.  Comparison  of  calculated  and  measured  results  demonstrate 
that  a  highly  uniform  optical  field  in  one  transverse  dimension  with 
efficient  utilization  of  optical  energy  can  be  obtained  by  using  the  OOF 
technique.  This  technique  can  be  used  for  pabticle  sizing  and  velocimetry. 


1.  Introduction 

A  uniform  intensity  distribution  across  a  laser  beam  is  important  in 
particle  measurements  using  Mie  theory.  Since  a  laser  beam  has  a  Gaussian 
intensity  distribution,  signals  scattered  by  a  small  particle  passing 
through  the  center  of  a  beam  and  by  a  larger  particle  passing  through  a 
region  of  lower  intensity  may  not  be  distinguishable.  Several  techniques 
have  been  developed  to  solve  this  problem,  which  can  be  divided  into  two 
groups.  In  the  first  group,  the  Gaussian  distribution  is  directly  used  to 
infer  properties  of  the  particle  size  and  velocity  distribution  [1-3],  in  the 
second  group,  the  Gaussian  distribution  is  first  transformed  into  some 
approximation  of  a  uniform  intensity  distribution  [4-6],  In  this  paper,  a 
novel  single-particle  measuring  technique  of  the  second  type,  designated 
as  the  Overlapping  Optical  Field  (OOF)  technique,  is  described.  This 
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Using  dynamic  holography  for  iron  fibers  with  submicron 
diameter  and  high  velocity 

6.W.  Xie,  L.F.  Hua,  S.  Patel,  P.  Scott,  D.T.  Shaw 

State  University  of  New  York  at  Buffalo,  Department  of  Electrical 
and  Computer  Engineering,  Buffalo,  New  York  M260 


ABSTRACT 

A  divergent  laser  beam  was  used  to  record  holograms  of  moving  iron  fibers  with  submicron 
diameter  ( 0.5  microns)  and  high  velocity  ( 1 25  meters  per  second).  Analysis  based  on  geometric  optics 
is  used  to  derive  the  magnification  factors  for  recording  and  reconstruction.  Experimental  results 
consisting  of  single-shot  holograms  of  fast  moving  fibers  are  presented.  This  technique  can  be  used  for 
dynamic  measurement  of  particle  size  and  fiber  orientation  distribution. 


1.  INTRODUCTION 

The  in-line  far-field  holographic  technique  has  been  used  for  particle  measurement  for  many  years, 
and  a  wide  variety  of  useful  results  have  been  obtained1  “ 1 8.  But  up  to  the  present,  there  is  no  report 
in  the  literature  demonstrating  dynamic  measurement  of  moving  fibers  particles  with  both  submicron 
diameter  and  high  velocity  (over  1  OOm/sec). 

A  difficulty  needing  to  be  overcome  in  this  case  is  that  dynamic  holography  of  very  small,  rapidly 
moving  particles  fields  requires  both  the  use  of  a  microscope  objective  lens  of  high  magnification  and 
long  working  distance,  and  high  power  pulsed  lassr.  For  in-line  holography,  a  plane  wave  is  the  obvious 
choice  to  illuminate  the  particle  field.  But  as  the  convergent  point  of  microscope  objective  lens  is 
always  close  to  one  element  of  the  objective  lens,  when  the  incident  wavefront  is  planar  the  extreme 
power  densities  casued  by  focusing  a  high  power  pulsed  laser  introduces  heating  of  the  lens  surface 
sufficient  to  damage  the  lens.  In  this  paper  we  will  present  a  divergent  beam  method  which  prevents 
damage  of  the  microscope  objective  lens.  The  pulse  energy  received  by  the  surface  of  the  last  piece  of 
the  microscope  objective  lens  is  spread  over  a  considerably  larger  crossection  than  the  planar  case. 
Analysis  by  geometric  optics  shows  that  an  approximately  6°  divergent  beam  spreads  the  energy  over 
60  times  as  large  a  surface  as  the  corresponding  planar  beam ,  shifting  the  incident  radiant  power 
density  from  well  above  the  amount  lethal  to  the  lens  to  safety  below. 

In  Section  2  the  divergent  beam  method  is  described  and  a  quantitative  analysis  of  the  power  density 
impinging  on  the  microscope  lens  over  the  range  of  beam  divergence  angles  from  0°  to  T  is  presented. 
In  Section  3  the  magnification  factors  for  recording  and  reconstruction  are  calculated.  Section  4 
describes  the  results  of  divergent  beam  in-line  holography  for  fast  moving  narrow  gauge  iron  fibers, 
and  includes  an  error  analysis.  A  summary  is  offered  in  Section  5. 

2.  DIVERGENT  LASER  BEAM  METHOD 

The  microscope  objective  lens  ( ML)  which  we  used  for  taking  dynamic  holograms  of  fine  fibers  is  a 
product  of  Bausch  &  Lamb  Corp.,  25x  0.3 1NA.  It  is  a  compound  objective  lens  which  consists  of  five 
lens  pieces:  ML] ,  ML2,  MLj,  ML4,  and  ML5.  Two  of  these  are  Fraunhofer  cemented  doublet  lens, 

ML  1  2,3,4  are  positive  lenses  and  ML5  is  a  negative  lens.  The  structure  of  ML  is  sketched  in  Fig.l, 


19 


Single  Shot  Holography  of  Fast  Moving  Fibers  with 
Submicron  diameter 

G.W.  Xie,  L.F.  Hua,  S.  Patel,  P.  Scott,  D.T.  Shaw 

State  University  of  New  York  at  Buffalo,  Department  of  Electrical 
and  Computer  Engineering,  Buffalo,  New  York  14260 


ABSTRACT 

A  divergent  laser  beam  was  used  to  record  holograms  of  moving  iron  fibers 
with  submicron  diameter  0.5|im  and  high  velocity  125  meter  per  second.  The 
principle  of  divergent  laser  beam  method  based  on  geometrical  optical  theory 
and  the  magnification  calculation  of  recording  and  reconstructing  process  of 
hologram  are  presented.  The  experimental  results  of  single-shot  holograms  of 
fast  moving  fibers  are  demonstrated.  This  technique  can  be  used  for  dynamic 
measurement  of  particle  size  and  fiber  orientation  distribution. 


1  INTRODUCTION 

The  in-line  far-field  holographic  technique  has  been  used  for  particle 
measurement  for  many  years  and  great  successful  achievements  have  been 
obtained  in  this  field1'18.  But  up  to  the  present,  there  is  no  report  in  the 
literature  about  the  dynamic  measurement  of  moving  fibers  particles  with  both 
submicron  diameters  and  high  velocity  (over  lOOm/sec).  The  reason  is  that 
dynamic  holograms  of  very  small  rapidly  moving  particles  fields  requires  both 
the  use  of  microscope  objective  lens  of  high  magnification  with  long  working 
distance  and  high  power  pulse  laser.  According  to  the  general  principle  of 
in-line  holographic  technique,  the  plane  wave  is  usually  used  to  record  the 
hologram.  But  as  the  convergent  point  of  microscope  objective  lens  is  always 
close  to  one  piece  of  the  lenses  which  compose  the  whole  objective  lens  when 
incident  is  plane  wave,  the  extreme  power  densities  derived  by  focusing  a  fast 
high  power  pulse  laser  introduces  heating  of  the  lens  surface  sufficient  to 
destroy  the  lens.  In  this  paper  we  will  present  a  divergent  beam  method  to  avoid 
the  damage  of  the  microscope  objective  lens,  By  this  method,  the  pulse  energy 
received  by  the  surface  of  the  last  piece  of  the  microscope  objective  lens  can  be 
spread  over  a  considerably  larger  crossection  than  the  planar  beam  is  used. 
Analysis  by  geometric  optics  shows  that  an  approximately  6°  divergent  beam 
covers  60  times  as  large  a  crossection  of  the  entire  lens  surface  as  the 
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toward  attending  this  Symposium  in  the  hope  of  meeting  you  and  hearing  of  your 
work.  If  you  suspect  that  you  are  unable  to  attend,  please  advise  us  immediately 
so  that  your  name  doesn't  appear  in  the  Advance  Program  or  in  the  Final  Program 
if  too  late. 


Should  you  need  any  information,  don't  hesitate  to  contact  Sylvie  Trahais  at 
ESIC  at  the  address  mentioned  at  the  top  of  this  letter.  Please  note  that  the 
paper(s)  number(s)  which  appears  on  the  above  label  is  to  be  recalled  in  all 
correspondance  (fax,  mail  or  even  telephone  calls...) 
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